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We consider the nearest-neighbor spacing distributions of mixed random matrix ensembles in- 
terpolating between different symmetry classes, or between integrable and non-integrable systems. 
We derive analytical formulas for the spacing distributions of 2 x 2 or 4 x 4 matrices and show 
numerically that they provide very good approximations for those of random matrices with large 
dimension. This generalizes the Wigner surmise, which is valid for pure ensembles that are recovered 
as limits of the mixed ensembles. We show how the coupling parameters of small and large matrices 
must be matched depending on the local eigenvalue density. 
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I. INTRODUCTION 

Random matrix theory (RMT) is a powerful mathe- 
matical tool which can be used to describe the statistical 
behavior of quantities arising in a wide variety of com- 
plex systems. It has been applied to many mathematical 
and physical problems with great success, see [iHS] for 
reviews. This wide range of applications is based on the 
fact that RMT describes universal quantities that do not 
depend on the detailed dynamical properties of a given 
system but rather are determined by global symmetries 
that are shared by all systems in a given symmetry class. 

In RMT the operator governing the behavior of the sys- 
tem, such as the Hamilton or Dirac operator, is replaced 
by a random matrix with suitable symmetries. One then 
studies statistical properties of the eigenvalue spectrum 
of such random matrices, typically in the limit of large 
matrix dimension. To compare different systems in the 
same symmetry class with RMT, the eigenvalues of the 
physical system as well as those of the random matri- 
ces need to be "unfolded" |4^. The purpose of such an 
unfolding procedure is to separate the average behavior 
of the spectral density (which is not universal) from the 
spectral fluctuations (which are universal). Unfolding is 
essentially a local rescaling of the eigenvalues, resulting 
in an unfolded spectrum with mean level spacing equal 
to unity. How the rescaling is to be done is not unique 
and may depend on the system under study. 

In this paper we focus on the so-called nearest-neighbor 
spacing distribution P(s), i.e., the probability density to 
find two adjacent (unfolded) eigenvalues at a distance 
s. This quantity probes the strength of the eigenvalue 
repulsion due to interactions and can be computed an- 
alytically for the classical RMT ensembles, resulting in 
rather complicated expressions given in terms of prolate 
spheroidal functions [5]. However, it was realized early 
on that the level spacing distribution of large random 
matrices is very well approximated by that of 2 x 2 ma- 
trices in the same symmetry classPj For most practical 



purposes it is sufficient to use this so-called Wigner sur- 
mise [HI p. 199] instead of the exact analytical result. It 
is given by 

Ppis)^aps^e-''^'' (1.1) 

with f3 — 1,2,4 corresponding to the Gaussian orthog- 
onal (GOE), unitary (GUE), and symplectic (GSE) en- 
semble of RMT, respectively. The quantities ap and &^ 
are chosen such that 

/•oo poo 

/ dsPp{s) = l and (s) = / ds Ppis) s 1 (1.2) 
Jq Jq 

in all three cases. Explicit formulas will be given in 
Sec.ini 

RMT describes quantum systems whose classical coun- 
terparts are chaotic and correctly predicts the strong 
short-range correlations of the eigenvalues due to inter- 
actions. In contrast, the level spacing distribution of a 
quantum system whose classical counterpart is integrable 
is given by that of a Poisson process, 

Pois)^e-', (1.3) 

corresponding to uncorrelated eigenvalues. We assign the 
Dyson index /3 = to ensembles of this kind, which is a 
consistent extension of the generalized Gaussian ensem- 
bles with arbitrary real f3 > introduced in [lO] . 

Often physical systems consist of parts with different 
symmetries, or of a classically integrable and a chaotic 
part. Changing a parameter of the system may then 
result in transitions between different symmetry classes. 
We assume the Hamiltonian describing the system to be 
of the forrr0 

H = Hfi + XHp, , (1.4) 

where -ff^ represents the original system whose symme- 
try/integrability is broken by the perturbation Hjsi for 



This does not work for non-Hermitian complex matrices (6] [7] . 



^ Other possibilities have also been investigated, see, e.g., [5l 1111 
I12| . but will not be considered in this paper. 
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small coupling parameter A, and vice versa for large A. 
For the quantities we analyze the absolute scale of H is 
irrelevant, only the relative scale between the different 
parts matters. 

From the level statistics point of view, Hp and Hj^i 
correspond either to a Poisson process or to one of the 
three RMT ensembles. Hence, there are (2) = 6 possi- 
bilities for a transition between two of these four cases 



in Eq. (1.4), i.e., Poisson-GOE, Poisson-GUE, Poisson- 



GSE, GOE-GUE, GOE-GSE, and GUE-GSE. If a GSE 
matrix is involved in the transition, there are two possi- 
bilities for the other matrix: self-dual or notl3 This leads 
to an even larger variety of mixed ensembles. Many tran- 
sitions of this kind have been studied in earlier works, 
usually for large matrix dimension. Transitions between 
Gaussian ensembles are considered in S], but closed 
forms for the spacing distribution could not be obtained, 
and self-dual symmetry was not conserved in the tran- 
sitions involving the GSE. Mixtures of Gaussian ensem- 
bles with conserved self-dual symmetry and small matrix 
size are considered in [13' , but only numerical results are 
given for the spacing distributions. Other examples in- 
clude the heuristic Brody distribution |14j interpolating 
between Poisson and the GOE, and a complete study of 
the transition between Poisson and the GUE [T5] . 

Note that an exact analytical calculation of P{s) for 
systems described by an Ansatz of the form ( 1.4 ) is much 



harder than, e.g., the analytical calculation of low-order 
spectral correlation functions, which are already difficult 
to obtain. Here, we do not attempt an analytical calcu- 
lation of P{s) for large matrix dimension. Rather, moti- 
vated by the reliability of the Wigner surmise, we study 
the possible transitions in Eq. ( 1.4 1 for 2 x 2 matrices (or. 



in the symplectic case, 4x4 matrices, because the small- 
est non-trivial self-dual matrix has this size) and compare 
the resulting level spacing distributions with that of large 
random matrices, the latter obtained numerically. The 
cases of Poisson-GOE and GOE-GUE were worked out 
earlier by Lenz and Haake [16] , and the spacing distribu- 
tion of a 2 X 2 matrix interpolating between Poisson and 
GUE is given in [T7] . These cases will briefly be reviewed 
below, and the remaining ones are the main subject of 
this work. 

This paper is organized as follows. In Sec. [H] we derive 
analytical results for P{s) for small matrix sizes. If Hj^i 
is from the GSE (i.e., Hpi is self-dual) we construct in 
Sees. |IID| |II F| and |II G| self-dual matrices Hp to main- 



tain the Kramers degeneracy. In Sec. II H we consider 
the case where a 4 x 4 GSE matrix is perturbed by a 
non-self-dual GUE matrix. Section |lTl] provides strong 
numerical evidence that the results obtained in Sec. |TT] 
approximate the spacing distributions of large random 
matrices very well. We give a perturbative argument for 



the matching of the couplings used for the Wigner sur- 
mise and for large matrices, respectively, and derive an 
approximate result that involves the eigenvalue density. 
This result describes the numerical data rather well. We 
also show that the transitions from the GSE to cither a 
non-self-dual Poissonian ensemble or the GOE proceed 
via an intermediate transition to the GUE and can also 
be described by the surmises calculated in Sec. |lT| We 



summarize our findings and conclude in Sec. IV Techni- 
cal details are worked out in several appendices. 



II. SPACING DISTRIBUTIONS FOR SMALL 
MATRICES 

A. Preliminaries 

In the spirit of the Wigner surmise, we now calculate 
the distributions P{s) of eigenvalue spacings s of mixed 
ensembles for the smallest nontrivial (i.e., 2 x 2 or 4 x 4) 
matrices, with P{s) normalized as in Eq. ( 1.2 ). Unfolding 
is not needed for these matrices since they have only two 
independent eigenvalues (except for Sec. II H). We first 



study the transitions from the integrable to the chaotic 
case for the three Gaussian ensembles and then proceed 
to the transitions between different symmetry classes. 
We define the 2x2 Poisson process by a matrix 



^0 = 




p 



(2.1) 



where p is a Poisson distributed non-negative random 
number with unit mean value, i.e., its probability den- 
sity is Pa{p) — e^P. The eigenvalue spacing of this ma- 
trix is obviously Poissonian, as the spacing is just p, and 
therefore we obtain Eq. (1.3). 



The choice of Hq may look like a special case, but it suf- 
fices for our purposes. The most general Hermitian 2x2 
matrix with spacing p can be obtained from Eq. (2.1) 



^ An even-dimensional matrix A is called self-dual if JA^J^ ■ 
with J given in Eq. ||F4[|. 



by a common shift of the eigenvalues (which does not 
influence the spacing) and a basis transformation. This 
transformation can be absorbed in the perturbing ma- 
trix since it does not change the probability distribution 
of the latter. To see this suppose we had started with a 
general nondiagonal TJqj also with eigenvalues and p, 
instead of Eq. (2.1 1. When added to a random matrix 
Hp' with /?' = 1, 2,4, we choose it to be real symmetric, 
Hermitian, or self-dual, respectively, in order to preserve 
the symmetry properties of Hp/. Then i?o is diagonal- 
ized by a suitable matrix $7, i.e., diag(0,p) — il^^Hofl, 
where is orthogonal (/3 = 1), unitary (/3 = 2), or sym- 
plectic (/3 = 4). In the total matrix H this is equivalent 
to il,~^Hj3>n perturbing diag(0,p), but the probability 
distribution of the perturbation is invariant under the 
transformation il. 

For matrices i?i,2,4 from the GOE, GUE, and GSE, 
respectively, we choose the mean values of the matrix 
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FIG. 1. Spacing distr ibutions Po^;9'(s;A) of the transitions Poisson GOE (left, Eq. (pl|), GUE (middle, Eq. < |2.12[ )) and 
GSE (right, Eq. (2.26l) for 2 x 2 or 4 x 4 matrices with coupling parameters A — 0.02, 0.08, 0.2, 0.4, and 0.8 (maxima moving 
from left to right as A increases). In the GSE case the matrix representing the Poisson process was made self-dual. All formulas 
were verified by comparison to numerically obtained spacing distributions of 2 x 2 or 4 x 4 random matrices. 



elements to be and the normalization 



1.2,4) 



[(-f^l!2,4A^jJ 



1, 



(2.2) 



The index i> ~ 0, . . . , (3 ~ 1 distinguishes the components 
of the complex/quaternion GUE/GSE matrix elements, 
while the GOE matrix elements possess only a real part. 



All results we derive from Eq. (1.4) will be symmetric 
in A since the distribution of the elements of H/3> is sym- 
metric about zero (the perturbation will be taken from 
one of the Gaussian ensembles in each case) . This means 
that our results should be expressed in terms of |A|. To 
avoid such cumbersome notation we restrict ourselves to 
non-negative A. 



B. Poisson to GOE 

We first consider the case that corresponds to a classi- 
cally integrable system perturbed by a chaotic part with 
anti-unitary symmetry squaring to 1. The integrable part 
is represented by a Poisson process, and the chaotic one 
by the GOE. The spacing distribution for this case has 
been derived in [T^], and we state it here for the sake of 
completeness. 

The 2x2 random matrix 



H = Ho + XHi = 




p 



+ A 



a c 
c b 



(2.3) 



consists of Hq from (2.1) and Hi from the GOE, i.e., 



a real symmetric matrix with normalization given in 
Eq. (2.2). The calculations are very similar to the ones 
for the transition from Poisson to the GSE, which are 
presented in Sec. |II D[ The resulting spacing distribution 
of H reads 



Fo^i(s;A) = Cse 



with 



DiX) 



2A 



dx ( 



xDs 



0,A2) 



(2.4) 



(2.5) 



C(A) = 2D{X)^ 



(2.6) 



where U is the Tricomi confluent hypergeometric func- 
tion (or Kummer function) [TBI Eq. (13.1.3)] and Iq is 
a modified Bessel function [181 Eq. (9.6.3)]. Po-i.i(s;A) 
is plotted in Fig. [l] (left) for various values of A. The 
formula is equivalent to the one given in [TB , but our 
integration variable x is scaled differently. 

In the limiting cases of A — > and A — > oo we have 



D{X) 



/l/(2A) 



for A 
for A 



0, 



(2.7) 



Using the asymptotic expansion of the Bessel function, it 
is straightforward to show that for A — > we obtain the 
Poisson result e"*. It is even simpler to show that the 
Wiener surmise (7rs/2)e-'^" for the GOE is obtained 
for A — oo. 

The small-s behavior of Pq^i{s; X) shows interesting 
features. To investigate this behavior, we consider sepa- 
rately the cases A = and A > 0. For A = we have by 
construction 



For A > we obtain from Eq. ( |2.4[ ) 

Po^i{s;X)^c{X)s + 0{s^) 

with 



c(A) 



2A 



for A ^ . 



(2.8) 
(2.9) 

(2.10) 



which means that we recover the linear level repulsion of 
the GOE for arbitrarily small A, i.e., for arbitrarily small 
admixture of the chaotic part as also observed in [T9H2T] . 
This implies that for A — >■ the distribution, viewed as 
a function of A, develops a discontinuity at s = 0, since 
Po^i(s = 0;A = 0) = 1 while Po^i{s = 0;A > 0) = 0. 
This effect is clearly seen in Fig. [l] (left). 

For small values of A and s, we observe something rem- 
iniscent of the Gibbs phenomenon, i.e., the interpolation 
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overshoots the Poisson curve considerably. In the Umit 
of A — 0, one can show (see App. |A 2[ ) that the maxi- 
mum of Po->i is at Smax — 2.51393 A with a finite value of 
Po->-i(smax; A — > 0) = 1.17516. This implies an overshoot 
of 17.5% compared to the Poisson curve. Such an effect 
also occurs in the transitions from Poisson to GUE and 
GSE that are treated in Sees. |II C] and [ITD| below, with a 
quadratic/quartic level repulsion in the small-s regime. 



matrix with normalization (2.2). The spacing distribu- 



tion of an equivalent setup with different normalizations 
of the random matrix elements was already considered in 
[T7] . so we just state the result, 



with 

Di\)-- 
C(A) = 



dxe 4^ 



sinh2^ 



(2.12) 



/= + -e^ erfc(A)-^Ei(A^) 



2A^ 



2F2 (5,1; |, |; A^) , 



4A 



D{Xf 



(2.13) 
(2.14) 



Here, erfc is the complementary error function [181 
Eq. (7.1.2)], Ei is the exponential integral [THl 
Eq. (5.1.2)], and 2F2 is a generahzed hypergeometric 
function [#1 Eq. (9.14.1)]. 

To check the validity of Eq. ( |2.12 ) and to see the emer- 
gence of the limiting spacing distributions, we now con- 
sider the limits A — t- and A ^ 00. First note that for 
A — we have 



D 



1 

2A 



and C ' 



1 



A^ 



(2.15) 



so that Eq. (2.121 becomes for s > 



Po^2(s;0) = lim 



dxe 4^ ^ ' 



sinh^ 



— -1= I dx lim — { e 

SyTT Jq X A^o A 



(2.16) 



which is the Poisson distribution as required. For A — > 00 
we have 

16A 

7^ 



D 



and C 



(2.17) 



so that Eq. (2.12 1 becomes 



Po^2(s; 00) = lim 

A— f 00 

32s2 



16Ase 



r3/2 



dx t 



sinh 4^ 



C. Poisson to GUE 

We now consider the transition from Poisson to the 
GUE. This corresponds to a classically integrable system 
with a chaotic perturbation without anti-unitary symme- 
try. The 2x2 random matrix 

//^h.+ah..(»»)+a(^^°,^_="V") 

contains H2 from the GUE, i.e., a complex Hermitian Eq. (|2.12 1 



(2.18) 



which is the Wigner surmise for the GUE. 

The integral in Eq. (2.12) can be computed numeri- 



cally without difficulties as the integrand decays like a, 
Gaussian for large x and becomes constant for small x 
The resulting distribution Po-i-2(s; A) is plotted in Fig. 
(middle). 

As in Sec. |IIB[ a discontinuity is found at s = to- 
wards the Poisson result. For A > we obtain from 



Po^2(s;A)=c(A)s2 + 0(s4) 



with 



c(A) 



1 
2A2 



for A ^ . 



(2.19) 



(2.20) 



Hence we obtain the quadratic level repulsion of the GUE 
for arbitrarily small coupling parameter. For A — > 0, the 
maximum of the function is at s„iax — 3.00395 A, with a 
value of Po^2(smax; A ^ 0) = 1.28475 (see App.[A2). 



D. Poisson to GSE 

In this classically integrable system is perturbed 

by a chaotic part with anti-unitary symmetry squaring 
to —1 and hence represented by the self-dual matrices 
of the GSE. One has to consider 4x4 matrices here, 
because a self-dual 2x2 matrix is proportional to I2 and 
has only one non-degenerate eigenvalue. As mentioned 
in the introduction, there are now two possibilities: The 
Poisson process could be represented by a self-dual or a 
non-self-dual matrix. Here we only consider the former 
possibility, while the latter will be discussed in Sec. |IIIEl 
A self-dual Poisson matrix is obtained by taking a tensor 
product of Ec^. (2.1) with I2. Thus the transition matrix 
is 



/O 0\ 



p 

\0 p/ 



* Note that the integral can be expressed in terms of imaginary 
error functions, but for increasing s dehcate canceUations occur 
that make it impractical to use this form for numerical evalu- 
ation. This is why we present Eq. | |2.12| as the final formula, 
which is well suited for numerical integration. 



5 



+ A 



a 




Co + ^C3 Ci + IC2^ 

— ci + ic2 Co — ics 



Co " «C3 -Ci - IC2 

\,ci — ic2 Co + ic3 



(2.21) 



where the GSE matrix is Hermitian and self-dual, and 
can be represented by a 2 x 2 matrix whose elements are 
real quaternions, see [5| for details. 

We now explain the calculation of the spacing distribu- 
tion for this transition. The computation of the previous 
cases, Poisson to GOE and Poisson to GUE, can be done 
in a similar fashion. 

Due to the self-dual structure of H, the spacing S be- 
tween its non-degenerate eigenvalues spacing can be com- 
puted analytically and reads 

5 = A^(a-6-p/A)2 + 4c^c^, (2.22) 

where the repeated index /i indicates a sum from to 
3. We have intentionally written S instead of s since we 
eventually need to rescale the spacing to ensure (s) = 1. 
The desired spacing distribution is proportional to the 
integral 

3 

I{S) - / dpdadb Y[dc, Po{p)Pa{a)Pb{b)Pc,Ac.) 

x5(^S-^[a-{b + p/A)]2 + 4c^c^^ , (2.23) 

where we have rescaled 5* by A for simplicity and are not 
yet concerned with the normalization. The distributions 
Pa{a) of the random variables a = a, b, cp , ci, C2, C3 are 
Gaussian, with variances given by Eq. (2.2), 



2 O/T^ 1 

a^h co,ci,C2,C3 



(2.24) 



Inserting this into Eq. (2.231 and shifting b ^ b — p/X 



gives 



I{S)(x dp dadbJldc^e'P-^''' 

Jo J —00 , n 



X S [ S - (a - b)^ + 'iCf.c^ 



(2.25) 



The multi-dimensional integral in this expression is com- 
puted in App. |B 1] Rescaling the spacing and normalizing 
the distribution to satisfy Eq. (1.2), we obtain 



Po^4(s;A) = Cs*e 

X / (1 - a;2) e("^")'+2^"^" erfc(a;Z?s -t- A) 



(2.26) 



with 



D{X) 



dx e 



-2Xx 



(2.27) 



(Ax^ + 2a:)e-^ + y/^{4x^ + 4x'^ - 1) erf (x) 



C(A) = 2Ae^'L»(A)' 



(2.28) 



where erf is the error function [THl Eq. (7.1.1)]. 
In the limiting cases of A — > and A 00 we find 



D{X) 



fl/(2A) forA^O, 
[8/(37^) for A ^00. 



(2.29) 



For A — > 0, the integral in Eq. (2.26) can be evaluated 



in saddle-point approximation (with the saddle point at 
X = —1), resulting in the Poisson result e~*. For A — )■ 00, 
we use the asymptotic expansion of erfc and obtain the 
Wigner sur mise (6 4/97r)3s4e-64^V9^ for the GSE. 

Equation (2.26) is plotted in Fig. [l] (right) and again 
displays a discontinuity at s = as A — >■ 0. For A > we 
now have 



with 



Po^4(s;A) =c(A)s* + 0(s<') (2.30) 
c(A)^^ forA^O. (2.31) 



For A — > 0, the maximum of the function is at Smax — 
3.76023 A, with a value of Po^4(smax; A ^ 0) = 1.43453 
(see App. |A2[). 



E. GOE to GUE 

With this subsection we start the investigation of tran- 
sitions between different chaotic ensembles using the 
smallest possible matrix size. 

We consider the 2x2 matrix 



H = Hi + XH2 



(2.32) 



The spacing distribution for this transition was already 
computed in I16| . With the normalization of ensembles 
given in Eq. (|2.2|), it reads 



Pi^2is;X)^Cse 



with 



DiX) 



V1 + A2 



A 

1 + A2 

C{X)^2Vl + X^D{Xf . 



arccot A 



(2.33) 

(2.34) 
(2.35) 



This formula matches the result of [H] up to a rescaling 
of the coupling parameter A by a factor of which is 
due to a different normalization of the ensembles used 
there. 

In the limiting cases of A — > and A — > cxi we have 



DiX) 



for A ^ , 
2/y^ for a — >■ 00 . 



(2.36) 
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FIG. 2. Spacing distributions Pg^g > (s; A) of the transitions GOE ^ GUE (left, Eq. (|2^), GOE ^ GSE (middle, Eq. ( [2ll| )) 
and GUE — >■ GSE (right, Eq. (2.50l) for 2 x 2 or 4 x 4 matrices with coupling parameters A = 0.1, 0.2, 0.4, and 0.8 (maxima 
increasing with A). In the cases involving the GSE, the GOE or GUE matrices were made self-dual. All formulas were verified 
by comparison to numerically obtained spacing distributions of 2 x 2 or 4 x 4 random matrices. 



For A — > 0, the error function in Eq. (2.33) can be re- 
placed by unity (for s > 0), and we obtain the Wigner 
surmise for the GOE. For A oo, using the first-order 
Taylor expansion of the error function yields the Wigner 
surmise for the GUE. 

The result ( |2.33[ ) is plotted in Fig. [2] (left). In the 
small-s region, we now have for A > 



l->2 



(s;A) = c(A)s2+0(s4) 



with 



c(A) 



2A 



for A ^ . 



(2.37) 



(2.38) 



Similar to the previous subsections, a non-analytic transi- 
tion between weaker and stronger level repulsion develops 
as A — 0, except that now there is no jump in the func- 
tion itself but rather in its derivative at s = 0. Therefore, 
the stronger level repulsion takes over immediately in the 
small s-regime, if A > 0. As we shall see below, this also 
happens in the remaining transitions, GOE to GSE and 
GUE to GSE, and seems to be a characteristic feature of 
the mixed ensembles. 



F. GOE to GSE 

As the GSE is involved in this transition, we need ma- 
trices of size 4x4. Again there are two possibilities: 
The GOE matrix could be made self-dual, or it could be 
non-self-dual (as it generically is) . Here we only consider 
the former case, while the latter case will be discussed in 
As in 



Sec. |IIIE 

by 



we define a modified GOE matrix 



/a c ON 

a c 

c 6 

\0 c 6y 



(2.39) 



with real parameters a, &, c. This matrix is self-dual, so 
we can add it to a matrix from the GSE without spoiling 



the symmetry properties of the latter. Thus we consider 
H = Hi®l2 + \Hi, (2.40) 



where Hi and i?4 are normalized according to Eq. \2.2\ . 
The eigenvalues of the sum are doubly degenerate and 
can be calculated easily due to self-duality. 



After some algebra (see App. B 2 ) we obtain for the 
spacing distribution of H 

Pi^4(s;A) = Cs4e-(i+2^')^'^' 

X / dx(l-a;2)e("^^)' [/o(z)-/i(z)], (2.41) 
Jo 

where z = (1 — x^)D^s^ ^ Jo and Ii are modified Bessel 
functions, and 



D{\) = 
C{X) = 



A-A^ + (l + A2)2arccotA 



29/2 



27rAVl + A2 
A2(l + A2)3/2i?(A)^ 



(2.42) 
(2.43) 



In the limiting cases of A — > and A — > oo we have 

0, 



y^/(23/2A) 
8/(3^/2^A) 



for A 
for A 



(2.44) 



For A — > 0, we use the asymptotic expansion of the Bessel 
functions to simplify the integral over x in Eq. (2.41 ) and 
obtain the Wigner surmise for the GOE. For A — oo, the 
exponential and the difference of the Bessel functions in 
the integral over x can be replaced by unity, and the 
Wigner surmise for the GSE follows trivially. 

The distribution Pi_j.4(s;A) is plotted for several val- 
ues of A in Fig. [2] (middle) and displays a continuous 
interpolation between the GOE and GSE curves. In the 
small-s region, the level repulsion is of fourth order for 
non-vanishing A. This is visible in the plots and can be 
shown by expanding Pi_j.4(s; A) for A > and small s. 



Pi^4(s;A) = c(A)s4 + 0(s6) 



(2.45) 
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with 



c(A) 



12A3 



for A . 



(2.46) 



G. GUE to GSE 



Again, due to the presence of the GSE, we have two 
possibihties for the GUE: self-dual or not. The former 
case is simpler and analyzed here, while the latter case 
will be considered in Sec. |II H[ We first have to clarify 
how to obtain a self-dual 4x4 matrix whose eigenvalues 
have the same probability distribution as those of a 2 x 2 
matrix from the GUE. In analogy to Sec. |IIF| one could 
try H2 I2, but the resulting matrix is not self-dual. 
Instead, we consider the matrix 



H2 




(2.47) 



with H2 given in Eq. ( |2.11[ ). The eigenvalues of are 
obviously equal to those of H2, but twofold degenerate. 
Interchanging the second and third row and column of 



H2, we obtain the matrix 



a 


Co - ici 





a 


Co -I- ici 




(2.48) 



which is self-dual and has the same eigenvalues as A 
matrix of this form was already introduced in [13j . 

The proper self-dual matrix for the GUE to GSE tran- 
sition is thus 



H = m 



XH4 



(2.49) 



with given in Eq. ( 2.21 1. The calculation of the corre- 



sponding spacing distribution proceeds in close analogy 
with the one presented in App. |B 2] and we find the closed 
expression 

P2^4(s;A) = Ce-(^^^)' 

X [2{Ds)^ - ^/^Dse'^^'^^ eifL{Ds)] (2.50) 

with the imaginary error function crfi(a;) = — zerf(ia;) 
and 



DiX) = 
C(A) = 



1 

2A3 



2 + A^ - A' 



arccsch A 



(2.51) 
(2.52) 



where arccsch is defined in [ISl Eq. (4.6.17)]. 

In the limiting cases of A — >■ and A — > 00 we have 



DiX) 



r2/(A^) 

\8/(3Ay^; 



for A 
for A 



0, 
00 . 



(2.53) 



For A — )■ 0, the asymptotic expansion of the second term 
in the square brackets of Eq. (2.50) yields —1. This can 



be neglected compared to the first term in the square 
brackets, which gives the Wigner surmise for the GUE. 
For A — ^ 00, Taylor expansion of the square brackets in 



Eq. (2.501 yields the Wigner surmise for the GSE. 

The resuh ( |2.50| ) is plotted in Fig. [2] (right). In the 
small-s region, we have for A 7^ 



i(s;A) =c(A)s4 + 0(s6) 



with 



c(A) 



256 
37r3A2 



for A ^ . 



(2.54) 



(2.55) 



H. GSE to GUE without self-dual symmetry 

In this section, we consider a matrix taken from the 
GSE whose Kramers degeneracy is lifted by a perturba- 
tion taken from the GUE without self-dual symmetry. As 
we shall see, this case also gives a surmise for other tran- 
sitions involving the GSE and another ensemble with- 
out self-dual symmetry. We will return to this point in 
SeciniEl 



1. General considerations 

The 4x4 transition matrix is 

H = Hi + XH2 (2.56) 

with Hi taken from the GSE and H2 from the GUE, 



both in standard normalization, Eq. (2.2). As H2 has no 



self-dual symmetry, the two-fold degeneracy of the GSE 
spectrum is removed and eigenvalue pairs split up. If 
the perturbation is small, there are two different spacing 
scales in this setup, as shown in Fig. |3] where the pertur- 
bation of two nearest-neighbor eigenvalues is sketched: 



H. 



~7K- 



~7K 



FIG. 3. Perturbation of GSE eigenvalues removing the de- 
generacy. 
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• ^i: The spacings between previously degenerate 
eigenvalues, which are of the same order of magni- 
tude as the coupling parameter for small couplings. 
They are formed by the two smallest/largest eigen- 
values of H . 

• The intermediate spacing, which is formed by 
the second and third largest eigenvalue of H . In 
the limit A — this is the original spacing of the 
GSE matrix Ha. 

The joint probability density of the eigenvalues of H 
is given, up to a rescaling, by [3 Eq. (14.2.7)] 

X [h{d2i)h{di3) + h{d32)h{dii) - hid3i)h{di2)] (2.57) 

with 

Aiei,02,03,ei)^Y[{0,-9,), (2.58) 



with 



i<j 

h{x) = xe-^'/^' 



(2.59) 
(2.60) 

(2.61) 



As we are only interested in spacings and thus in differ- 
ences of eigenvalues, we introduce new variables 



tl — d21 =02 — 01: 
h = ^32 = 03 — 02 : 



(2.62) 
(2.63) 
(2.64) 



and keep the original variable 0i. The Jacobi determi- 
nant of this transformation is 1, and we can now perform 
the 01 integration, which results (up to a constant factor) 
in 



X exp 



:A(-ti, 0,^2,^2+ is) 
{ti+2t2 + hf + 2t\ + 2tl 



(2.65) 



X [h{ti)h{t3) - h{ti+t2)h{t2+t3) + h{ti+t2+t3)h{t2)] . 

We now derive the distributions of the two different kinds 
of spacings from this formula. We assume 0i < 02 ^ 
03 ^ 0i and include the resulting combinatorial factor of 
4! explicitly. 



2. Spacings between originally degenerate eigenvalues 

To obtain the distribution of the spacing between the 
two smallest eigenvalues of H (the two largest ones give 
the same result due to symmetry), we set ti = Si and 
integrate over t2 and ^3 from to 00. This results in the 
spacing distribution 

POO 

Pl^2{si\\) = CD dt2dt3P{DsiMM) (2.66) 
"'0 



C(A) = I ^-3/2 ;^-6 (2 + A2)', 



(2.67) 



D{\) = C{\) dSidt2dt3SiP{SiMM)- (2.68) 
Jo 

We replaced Si by si to indicate that this is the spacing 
on the unfolded scale, i.e., with a mean value of 1. One of 
the integrals could in principle be done analytically, but 
this results in such a lengthy expression that it seems 
more sensible to evaluate all integrals numerically. 

The distribution in the limit A — ?> can either be ob- 
tained by perturbation theory, see App. |C 1[ or by di- 
rectly evaluating the spacing distribution in the limit 
A 0. First note that 



lim 

A->-0 



rA3 



X h{x) — 6{x) , 



(2.69) 



where the A-dependence of /i, which is suppressed in our 
notation, plays a crucial role. As the mean value of the 
spacing Si on the original scale has to become arbitrarily 
small in the GSE-limit due to the Kramers degeneracy, 
we consider a rescaled spacing si = Si/X. Therefore 
h{Si) becomes for small A 



h{Si) = h{Xsi) « Asi e 



(2.70) 



With these considerations we obtain from Eq. (2.65) 



P{XSi,t2,t3) oc Ae-i[(^^^i+2*^+*^')'+2^'^^?+2*3] (2.71) 



2s2 
^/^A 



'1 S{t3){XSi + t2)iXSi +t2+ h)t2{t2 + ta) 



- ^(ASi + ta) ^(^2 + t3) Xsit2t3{Xsi +t2 + ts) 
+ 6{XSi+t2+t3)6{t2)XSi{XSi+t2){t2+t3)t3 



as A — >■ 0. The last two terms in square brackets vanish 
upon evaluation of the t2 and integrals, because the 
zeros of the arguments of their (5-functions lie outside of 
the integration region. Performing the ^3 integration in 
the first term we obtain for nonzero A and si 



A^2(Si;A) 



.2 -s? 

oc Sn e 1 



(2.72) 



Up to normalization and rescaling this is the spacing dis- 
tribution of a 2 X 2 GUE matrix. 



In the opposite hmit A — 7> 00 the result (2.661 reduces 



to the distribution of the first and last spacings of a pure 
4x4 GUE matrix. This distribution can be obtained from 
similar considerations, starting from [5j Eq. (3.3.7)]. 

The result ( 2.66 1 is shown in Fig. [4] (left and middle) for 
several values of A, along with the limiting distributions 
for A — >■ and A — >■ cxd. All these curves are very similar 
and can only be distinguished by the naked eye in the 
zoomed-in plot. 



We have validated the result (2.661 by comparing it 



to the spacing distribution of numerically obtained 4x4 
random matrices. 
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1 

0.8 
0.6 
0.4 
0.2 
0, 



-2x2 GUE 
4x4GUE(s^) 
-interpolations 



0.5 



1.5 



2.5 



0.94 



0.93 



0.92 



0.91 



0.9 









-- -2x2 GUE 

-4x4 GUE (s^) 
— interpolations 





0.9 



0.95 



-GSE 

- - 4x4 GUE (Sg) 

— Interpolations 




FIG. 4. Spacing distributions for the transition from GSE — > GUE (non-self-dual) for 4x4 matrices and various values of 
the coupling parameter A. Left: spacings si between previously degenerate eigenvalues. Middle: also si, but zoomed in to 
the rectangular region indicated in the plot on the left; 2x2 GUE stands for A — >■ 0; 4 x 4 GUE stands for A — > <x; for the 
interpolation curves we chose A = 0.4, 1, 2 (maxima decreasing). Right: intermediate spacing S2; GSE stands for A = 0; 4 x 4 
GUE stands for A — > 00; for the interpolation curves we chose A = 0.05, 0.15, 0.3, 1 (maxima decreasing). 



Perturbed GSE-spacing 



III. APPLICATION TO LARGE SPECTRA 



We now consider the perturbed spacing of the original 
GSE matrix, which was formed by the two degenerate 
eigenvalue pairs of H^. The distribution of this spacing 
is obtai ned b y setting t2 = S2 and integrating P defined 
in Eq. (2.651 over ti and ^3 from to 00. With proper 



normalization as given in Eq. (2.2), this yields 



CD dhdhP(h,Ds2M) (2.73) 



In this section we will show numerically that the for- 
mulas derived in Sec. HI] for small matrices describe the 
spacing distributions of large random matrices very well. 
This observation should be viewed as our main result. 

When comparing the results obtained from large ma- 
trices to our generalized Wigner surmises, a natural ques- 
tion is how the corresponding coupling parameters, i.e., 
A in Eq. (1.4 1, should be matched. This question will be 



with 



addressed in the next subsection based on perturbation 
theory, while the numerical results will be presented in 
the remaining subsections. 



C(A) = ^7r-3/2 ;^-6 (2 + A2)' 



(2.74) 



D{\) = C{\) dS2dtidhS2P{ti,S2,t3). (2.75) 
Jo 

Again, the replacement of ^2 by S2 means that this is the 
intermediate spacing on the unfolded scale, i.e., with a 
mean value of 1. 
In the limit A 



reduces to the 
le opposite limit 



the result ( |2.73 
Wigner surmise for the GSE, while in t 
A — ^ CX3 it reduces to the spacing distribution of the in- 
termediate spacing of a pure 4x4 GUE matrix, which 
can again be obtained from similar considerations. 

The result (2.73) is shown in Fig. [4] (right) for sev- 
eral values of A, along with the limiting distributions for 
A — > and A — 00. The maximum of the interpolation 
first drops down as A is increased from 0, while at a value 
of A around 1 it starts to rise again as the distribution 
approaches its A — >■ 00 limit. Note that the limiting dis- 
tributions of si and S2 for A — > 00, i.e., the red dashed 
curves in Fig. |4j turn out to be almost identical to each 
other and to the Wigner surmise for the GUE. 



We have also validated the result (2.73) by comparing 
it to the spacing distribution of numerically obtained 4x4 
random matrices. 



A. Matching of the coupling parameters 

The setup is most easily explained by means of the 
transition from Poisson to the GUE. The Poisson case 
is represented by a diagonal N x N matrix Hq with in- 
dependent entries 6i (i = 1, . . . , N), each distributed ac- 
cording to the same distribution V{0), which we choose 
independent of N. The eigenvalue density of Hq is thus 
Pa{0) — NViO), and the local mean level spacing is 
l/po{6). We consider 



OlH-) 



(3.1) 



where H2 \s a.n N x N random matrix taken from the 



GUE, subject to the usual normalization, Eq. (2.2). 

As in the 2x2 case, the eigenvalues 9i will experience a 
repulsion through H2 ■ We will show in first-order pertur- 
bation theory that the relevant quantity for the repulsion 
is a combination of the eigenvalue density of Hq and the 
variance of the matrix elements of i?2. 

Ordinary perturbation theory in a yields a first-order 
eigenvalue shift of the 9i of 



(1) 



a(i?2)» 



(3.2) 
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This shift does not lead to a correlation of the eigen- 
values, as it just adds an independent Gaussian ran- 
dom number to each of them. Therefore, the eigenval- 
ues remain uncorrelated, and their spacing distribution 
remains Poissonian. 

However, if there is a small spacing of order a between 
twc|^ adjacent eigenvalues 6k and 9g of Hq, first-order 
almost-degenerate perturbation theory |23j predicts that 
the perturbed eigenvalues are the eigenvalues of the ma- 
trix 



Ok 

9e, 



a 



{H2)kk 

{H2)ek 



{H2)ki 

(H2)u 



(3.3) 



This matrix is almost identical to the 2x2 matrix con- 
sidered in Sec. IIC Eq. (2.11), with two differences: (i) 
The unperturbed eigenvalues 6k and 6i are shifted, but 
this does not affect the spacing distribution, (ii) The 
mean spacing of the unperturbed eigenvalues is not 1, 
but 1/ pq{0). We dropped the subscript on the eigenvalue 
6 here, because adjacent eigenvalues are very close for 
large iV, and therefore po{6k) ~ Poi^e) — Po{&)- 

To be able to match to the 2x2 formulas, we have to 
correct for the different mean spacing of the unperturbed 
matrix. We can do this by multiplying the matrix in 



Eq. (3.31 hy pq{9) without affecting the normalized spac- 



ing distribution. This results in the relation 
X{6)^Po{e)a 



(3.4) 



between the coupling parameters of the 2x2 and the 
NxN case. Note that the 2x2 parameter A has acquired 
a dependence on the eigenvalue 9 oi H through the local 
eigenvalue density of Hq. To be able to describe the spac- 
ing distribution of H in the spectral region around 6 by 
the generalized Wigner surmise, we assume that we have 
to insert this X{d) into the 2x2 formulas. This choice 
of universal coupling parameter is in line with an "un- 
folded" coupling parameter mentioned in [151 125 s-'^d a 
similar result from perturbation theory |25) . Appendix [d| 
contains a calculation for large matrices in second-order 
perturbation theory, also showing that the strength of 
the perturbation to be used in the generalized Wigner 
surmise only depends on the combination po{9) a. 

We now turn from the example "Poisson to GUE" to 
the general case, which we write as 



H = H, 



aH. 



(3.5) 



The same considerations hold with two modifications: (i) 
The unperturbed matrix is not necessarily diagonal by 
construction. However, it can be diagonalized by a trans- 
formation that can be absorbed in the perturbation]^ We 



can therefore treat it as diagonal (with eigenvalues corre- 
lated as dictated by the unperturbed ensemble), (ii) The 
mean spacing of the unperturbed 2x2 (4x4) matrix 
from Sec. is 



So = 1 (Poisson) , si -^/tt (GOE) , 



S2 = ^ (GUE) , Si 



16 



(GSE) 



(3.6) 



Therefore, we now have to multiply Eq. ( |3.5[ ) by si3pp{9) 
to get the correct mean spacing for the unperturbed 
matrix. This results in a universal, but 0-dependent, 
coupling parameter 



X{6) = sppp{9)a 



(3.7) 



with the eigenv alue density Pi3(9) of the unperturbed ma- 
trix. Equation (3.7) holds for all the transitions we con- 



sider, and in each case P is the Dyson index of the un- 
perturbed ensemble. 

In turn, this perturbative argument provides us with a 
formula of how to choose the coupling a in large matrices 
in order to approximate the spacing distribution of H by 
2 X 2 (4 X 4) formulas with parameter A, i.e.. 



A 



(3.8) 



where Pi3{6) is the eigenvalue density in the spectral re- 
gion we wish to study. In this way we can choose a value 
of A resulting in a spacing distribution roughly in the 
middle of the two limiting cases. Choosing a in Eq. (3.5) 



without this guidance is likely to result in a spacing dis- 
tribution that is dominated by one of the limiting cases. 



B. Transitions from integrable to chaotic 

1. Check of Wigner surmise 

We first consider transitions from Poisson to RMT for 
matrices with N non-degenerate eigenvalues. The ex- 
plicit numerical realization is the Hamiltonian 



H^Ho 



A 



Po(0) 



(3.9) 



where H/^i is a matrix taken from one of the Gaussian en- 
sembles, with normalization as given in Eq. (2.2). Hq is 
the same matrix as in Eq. (3.1 1 for the perturbation Hpi 



in GOE or GUE, whereas a self-dual Hq is constructed 



For small a, we are unlikely to find three or more small (i.e., of 
order a) consecutive spacings. 
® Note that we choose the perturbations H^/ such that their prob- 



ability distribution is always invariant under the transformations 
that diagonalize Hfj, just like in the Poisson to RMT cases. How- 
ever, this does not work for some of the transitions between the 
GSE and ensembles without self-dual symmetry, which we dis- 
cuss separately in Sec. IIII El 
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FIG. 5. Spacing distributions for the transition of large matrices from Poisson to GOE (top), GUE (middle), and GSE 
(bottom) with several values of the coupling A in Eq. ( |3.9[ ). The histograms show the numerical data, while the full curves are 
the analytical results for 2 x 2 (4 x 4) matrices with fitted coupling parameter A, see Sees. |II B through II D The quantity A2 
defined in App.[E]is a measure of the fit quality, which is small for a good fit. Each plot has been obtained by diagonalizing 
50, 000 matrices with 400 non-degenerate eigenvalues. 



by a direct product with I2 as in Sec. |II D| if the per- 
turbation is taken from the GSE. We choose a Gaus- 
sian for the distribution of the eigenvalues of Hq, i.e., 
V{e) = (l/y27r)exp(-6'V2), so po(0) = N/V2n. From 



Eq. (3.8) we would then expect the spacing distribution 
in the center of the spectrum around = to be approx- 



imat ed by the corresponding 2x2 formulas (2.4), (2.12), 
and (2.26) with coupling A = A. 



As can be seen in Fig. [5] the formulas for the 2x2 ma- 
trices indeed describe the spectra of large matrices quite 
well in a wide range of the coupling parameter A. The 
spacing distribution was evaluated in the center of the 
spectrum, defined as the interval (—0.2,0.2), because the 
eigenvalue density is almost constant and equal to po(0) 
in this region so that no unfolding is needed. The analyt- 
ical curve was obtained by a fit (see App. |e] for details) 
of the 2x2 (or 4x4) formula to the numerical data 
with fit parameter A. As expected by the perturbative 
considerations, A comes out on the same order of mag- 
nitude as A, and almost matches for small A. However, 
A is considerably smaller than A for stronger couplings. 
Presumably, the repulsion of the many other eigenvalues 



in the spectrum not present in the smallest matrices has 
a squeezing effect on the spacing, which works against 
the repulsion caused by the perturbation. This would 
explain the smaller coupling parameter. 



2. Dependence of coupling parameter on eigenvalue density 



The considerations in Sec. |III A] imply a linear relation 
between local eigenvalue density and effective coupling 
parameter, Eq. (3.7), for matrices of the form given in 
Eq. (3.1). This means that a perturbation should have 



a different impact on the spacing distribution of a sin- 
gle matrix in different regions of its spectrum (as qual- 
itatively observed in [H]). This subsection provides a 
detailed analysis of this phenomenon. 

Again, we consider a diagonal Poissonian matrix Hq of 
large dimension perturbed by a matrix taken from one of 
the Gaussian ensembles Hp>, 

H ^ Ho + aHfj, . (3.10) 

This time we will choose some fixed a and look separately 
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FIG. 6. Transitions Poisson — >■ GOE (left), GUE (middle), GSE (right). Top: unperturbed and perturbed eigenvalue density, 
the latter obtained numerically. Bottom: efTective coupling A obtained from fits of the 2 x 2 (4 x 4) level spacings Pq^^i{s) 
(see Eqs. (2.41, (2.121, and (2.26l) as a function of the local eigenvalue density in 35 equally large windows of the spectrum. 
Linear fit and proportionality factor with errors defined by the 95%-confidence interval are given in the plots. The quantity 82 
defined in Eq. (3.12 1 is a measure of the fit quality, which is small for a good fit. The numerical data were obtained from 10^ 
random matrices of dimension 600 (GOE, GUE) or 800 (GSE). 



at difFerent parts of the spectrum of H with a varying 
eigenvalue density. According to Eq. ( |3.7[ ) the effective 
2x2 (4x4) coupling parameter A should be the product 
of a and the local eigenvalue density of Hq. In App. [D] 
we show in perturbation theory up to second order that 
the local coupling parameter is in fact a function of this 
product. 

To treat such a system numerically one has to con- 
struct a Poissonian ensemble with a varying eigenvalue 
density, perturb it, and measure the coupling parameter 
in difFerent parts of the spectrum. This is done by cut- 
ting the spectrum into small windows with approximately 
constant eigenvalue density and fitting (see App. |e] for 
details) the spacing distributions inside the windows to 
the formulas for the 2x2 (4x4) matrices. We therefore 
obtain a fitted coupling parameter A for each window. 

For the numerical calculations, the eigenvalues 9i of the 
matrix Hq were distributed in the interval {—N/2, N/2) 
according to the somewhat arbitrarily chosen distribution 



-Pie,) 



1 

N 



N 



N 



(3.11) 



N being the number of independent eigenvalues of Hq 
The matrix Hpi is normalized in the usual way, Eq. (2.2 1 



The eigenvalue density p{9) of the total matrix H is 
plotted along with the analytical po{6) = N'P{9) of Hq in 
the top row of Fig. [6] One can see that the perturbation 



only has a negligible effect on the spectral density. 

The dependence of the coupling parameter on the 
eigenvalue density is plotted in the bottom row of Fig. [6] 
for a = 0.1. No error bars are shown because the statis- 
tical errors are negligibly small. A linear fit through the 
origin with minimized squared deviation was performed 
to obtain the proportionality factor between the eigen- 
value density and the coupling parameter. The quantity 
S2 shown in the plots is a measure of the fit quality and 
defined by 



\ 



N 

E 



(A. - 

N 



^ A- 



N 



(3.12) 



where the A^ are the numerically obtained coupling pa- 
rameters for each spectral window and the Ai are the 
corresponding predictions from the linear fit at the given 
eigenvalue density. Because 62 is a monotonically increas- 
ing function of the squared deviation it is also minimized 
by our fitting procedure. 

As can be seen, the linear dependence of the effec- 
tive coupling parameter on the eigenvalue density is con- 
firmed very well by the numerical data for all the transi- 
tions. Note that the fit quality gets better with increasing 
Dyson index /?', i.e., it is worst for the GOE and best for 
the GSE. This is most likely explained by the fact that 
the spacing distributions change more rapidly with re- 
spect to the coupling parameter for larger /3' (cf. Fig. [I]) , 
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FIG. 7. Spacing distributions for the transition of large matrices: GOE — > GUE (top), GOE — ^ GSE (middle), and GUE — >• 
GSE (bottom), with several values of the coupling A in Eq. (3.131. The histograms show the numeri cal d ata, whi le th e full 
curves are the analytical results for 2 x 2 (4 x 4) matrices with fitted coupling parameter A, see Sees. HE through II G The 
quantity A2 defined in App. |E] is a measure of the fit quality, which is small for a good fit. Each plot has been obtained by 
diagonalizing 50, 000 matrices with 400 non-degenerate eigenvalues. 



which allows for a more precise measurement of the cou- 
pling. 

Although the linear dependence of the effective cou- 
pling on the eigenvalue density has been demonstrated 
beyond reasonable doubt, the proportionality factor is 
less clear. As can be read off from Fig. |6] the proportion- 
ality factor is smaller than a, i.e., the measured coupling 
parameter is smaller than the expected one. This agrees 
with the observation in Sec. |III BTj where an explanation 
was given in terms of the effect of other eigenvalues. 



C. Transitions from one symmetry class to another 

1. Check of Wigner surmise 

We now consider chaotic systems composed of different 
symmetry classes, the latter represented by pure Gaus- 
sian ensembles. If the GSE is involved, we consider the 
case of a self-dual perturbed ensemble in this section (see 
Sec. |IIID] for the case of a non-self-dual perturbed ensem- 
ble). A self-dual GOE can be constructed by taking the 



direct product with I2 as in Se c. |IIF[ while the self-dual 
GUE is more involved, see App.|F[ All ensembles are nor- 
malized as in Eq. (2.2). Again motivated by Eq. (3.8 1, 



the Hamiltonian under consideration is 



H = Ha 



A 



-Ha 



(3.13) 



For large matrix size, the eigenvalue density of Hp is 
a semicircle which extends to = ^/2/3N, and its eigen- 
value density in the center is 



(3.14) 



The results for the three transitions among the Gaus- 
sian ensembles are shown in Fig. [7] for N = 400. Again, 
only the center of the spectrum, defined as the inter- 
val (—5,5), was evaluated (the whole semicircle extends 
to about ±28 for Hp e GOE and about ±40 for Hp G 
GUE). The coupling parameter A was obtained by a fit 
(see App. [e| for details) to the corresponding 2x2 (4x4) 
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formula, which yields a good approximation to the nu- 
merical data throughout the transition in each case. As 
in Sec. |IIIB 1[ A is close to A as expected. 

For A = 1 and N = 400 the mixed matrix is roughly 
given hy H — Hp + O(10^^)i?^'. From the A- values 
given in Fig. [7| which should be compared to those in 
Fig. [2j we see that the transition is almost completed 
in this case and that the spacing distribution is already 
very similar to the one of the perturbing ensemble. What 
is relevant for the transition is not the relative average 
magnitude of the matrix elements (which depends on N 
through the local eigenvalue density) but the rescaled 
coupling parameter A, i.e., the transition occurs for A = 
0(1). The same phenomenon was found for the two-point 
function '24' , which is related to the spacing distribution 
for small s. 



2. Dependence of coupling parameter on eigenvalue density 

We now consider the dependence of the coupling pa- 
rameter on the local eigenvalue density as in Sec. |III B2\ 
but now for transitions between Gaussian ensembles. In 
these cases, the fitting procedure of the effective coupling 
becomes less precise, because the functions of the spac- 
ing distributions change only very slowly with A, as can 
be seen in Fig. [2] Therefore, we restrict ourselves to the 
case of a self-dual GOE matrix Hi that is perturbed by 
a GSE matrix H4 as the level repulsion differs the most 
in these two ensembles. 

In Fig. [8] we show results from the mixed matrix 



H = Hi+ aHi , 



7 



Pi(0)si 



(3.15) 



with 7 — 0.2 (for details about 
and the normalization, see Sec 

Eq, 

be aX^si = 7p(0)/pi(O), 



the 



self-dual GOE 
IIICll). According to 



(3.8|) the effective 4x4 coupling parameter A should 

proportional to the lo- 



i.e., 

cal eigenvalue density normalized by the density /9i(0) in 
the center, with proportionality factor given by the input 
parameter 7. As one can see, there is again a linear de- 
pendence of the fitted coupling parameter A on the local 
density (and again, the perturbation has no measurable 
effect on the eigenvalue density). The proportionality 
factor is almost compatible with the expected value 7. 



D. Perturbation of a GSE matrix by a 
non-self-dual GUE matrix 

In this section, we apply the formulas derived in 



Sec. II H for the spacing distributions of a 4 x 4 matrix 



from the GSE perturbed by a matrix from the GUE, this 
time without self-dual symmetry, to large matrices. We 
consider a 2iV x 2A^ matrix 



H ^ Ha 



A 



P4(0)S4 



H. 



(3.16) 
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FIG. 8. Transition GOE ->• GSE. Top: unperturbed eigen- 
value density (approximated by a semicircle) and perturbed 
eigenvalue density. Bottom: effective coupling A obtained 
from fits of the 4x4 level spacings Pi_>4(s), see Eq. ( |2.41[ |, 
as a function of the local eigenvalue density in 35 equally 
large windows in the spectrum. Linear fit and proportional- 
ity factor with errors defined by the 95%-confidence interval 
are given in the plots. The quantity ^2 defined in Eq. (3.121 



is a measure of the fit quality, which is small for a good fit. 
The numerical data were obtained from 10^ random matrices 
of dimension 800. 



where H^ is taken from the GSE and H2 is the perturba- 
tion from the GUE. Both H^ and H2 are normalized in 



the usual way, see Eq. (2.2 1, and for the pref actor of H2 
see Sec. Ill A To ensure a constant eigenvalue density, 
we again restrict the measurements to the center of the 
spectrum, defined by the interval (—5,5). The numer- 
ically obtained spacing distributions were rescaled to a 
mean value of 1. 
As in Sec. 



II H 



we will separately consider the spac- 
ings between originally degenerate eigenvalues and the 
remaining ones. The distributions of the former were ob- 
tained by measuring every second spacing, starting with 
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FIG. 9. Spacing distributions between previously degenerate eigenvalues si (top) and previously non-degenerate eigenvalues 
S2 (bottom) for the transition GSE — >■ GUE without self-dual symmetry for various values of the coupling parameter A in 
Eq. (3.161. The histogra ms sho w the numerical data, while the full curves are the 2x2 GUE surmise P2 (top) and the surmise 
Pi-^2{s2', X) given in Eq. ( 2.73[ l (bottom), the latter with fitted coupling parameter A. The quantity A2 defined in App. [e| is 
a measure of the fit quality, which is small for a good fit. The numerical data were obtained by diagonalizing 50, 000 random 
matrices of dimension 400 for each plot. 



the first one of each random matrix. They are plotted in 
Fig. [9] (top) and show perfect agreement with the 2x2 
GUE surmise, which is practically indistinguishable from 
the exact result derived for 4x4 matrices, Pl^2{^i' ^) 
given in Eq. (2.661. As in Sec. II H 2 this distribution is 



almost independent of the coupling parameter. 

The distribution of the spacings between previously 
non-degenerate eigenvalues is shown in Fig. [9] (bottom). 
Again, every second spacing was measured, but starting 
with the second one this time. We get an almost per- 
fect agreement of the numerical data with the surmise 
P|^^, (s2;A) defined in Eq. (2.731 throughout the tran- 
sition. The parameter A was again determined by a fit 
(see App. [E]) and approximately matches the perturba- 
tive prediction from Sec. |III A| 



E. Other transitions between the GSE and 
ensembles without self-dual symmetry 

Let us now consider the transition from the GSE to 
either the GOE or Poisson, both without self-dual sym- 
metry. These two cases are more complicated than the 
cases discussed so far because, as we shall discuss now, 
the transitions proceed via an intermediate transition to 
the GUE. 

Let us first focus on the case 



H = Ha 



A 



P4(0)S4 



Hi, 



(3.17) 



where is from the GSE, Hi is from the GOE with- 
out self-dual symmetry, and we again concentrate on the 
central part of the spectrum (near zero). For small A, we 
show in App. [C] in first-order perturbation theory that 
the perturbation by the GOE has exactly the same effect 
on the eigenvalues as the perturbation by the GUE con- 
sidered in Sec. |IIII[ modulo a rescaling of the coupling 



parameter, i.e., 

PUi{si;\) = ^4^2(^1; A/\/2) ~ P2(si) 
PLi{s2;X)^PL2is2;X/V2). 



(3.18) 
(3.19) 



Therefore, we first expect a transition from the GSE to 
the GUE, corresponding to the breaking of the self-dual 
symmetry. This expectation is confirmed in Fig. [lO] (top 
and middle). 

As A is increased to very large values, a transition to 
GOE behavior must eventually occur. The question is 
whether this transition is described by the surmise of 
Sec. HE We show in Fig. 10 (bottom) that this is indeed 
the case. Note that a rising A amounts to a shrinking 
fitted coupling parameter A because the direction of the 
transition is turned around compared to Sec. |IIE[ Here, 
A — !■ cx) means that _ff is a pure GOE matrix, which is 
described by the surmise with A = 0. 

The case of GSE to Poisson without self-dual symme- 
try is analogous. For small values of the coupling param- 
eter, the self-dual symmetry of the GSE is broken by the 
perturbation so that we expect a GSE to GUE transi- 
tion for the spacings si and S2 as in the GSE to GOE 
case considered above. For very large values of the cou- 
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FIG. 10. Spacing distributions for the transition GSE — >■ GOE without self-dual symmetry for various values of the coupling 
parameter A in Eq. 13.171. Top: spacings si between previously degenerate eigenvalues (for small A). Middle: spacings S2 
between previously non-degenerate eigenvalues (also for small A). Bottom: all spacings (for large A). The histograms show the 
numerical data, while the full curves are the 2x2 GUE surmise P2 (top), the surmise P4_n(s2; A) given in Eq. ( |3.19 1 (middle), 
and the surmise Pi_>2(s; A) given in Eq. ( |2.33[ ), the latter two with fitted coupling parameter A. The quantity A2 defined in 
App. |E]is a measure of the fit quality, which is small for a good fit. The numerical data were obtained by diagonalizing 50, 000 
random matrices of dimension 400 for each plot. 



pling parameter we should eventually find a transition to 
Poisson behavior, described by the surmise of Sec. |II C| 
We have confirmed these expectations numericaUy but 
do not show the corresponding plots here. 

Note that in the transitions considered in Sees. IIIIBI 



through III D a single anti-unitary symmetry (or integra- 
bility in the case of Poisson) was broken or restored. In 
contrast, we now have two transitions. As A increases 
from zero, an anti-unitary symmetry T with — —1 



GSE 
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broken 



GUE 



rp2 
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A 







A~ 1 



FIG. 11. Schematic picture of the transition from GSE to non- 
self-dual GOE, which proceeds via an intermediate transition 
to the GUE. An analogous picture applies to the transition 
from GSE to non-self-dual Poisson. 



gets broken. As A decreases from infinity, either an 
anti-unitary symmetry with = 1 gets broken (in the 
case of GOE) or integrability gets broken (in the case of 
Poisson). For intermediate values of A the system fol- 
lows GUE statistics because all anti-unitary symmetries 
and/or integrability are broken. This is illustrated in 
Fig. [11] 



IV. SUMMARY 

We have derived generalized Wigner surmises for the 
nearest-neighbor spacing distributions of various mixed 
RMT ensembles from 2x2 and 4x4 matrices. If the 
GSE was involved in the transition, we have distinguished 
two cases: (i) perturbations of the GSE by a self-dual 
ensemble, and (ii) perturbations of the GSE by a non- 
self-dual ensemble, for which we separately considered 
two different kinds of spacings. 

We have shown that all of these distributions yield a 
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good description of the spectra of large mixed matrices 
when restricted to a range of constant spectral density. 
The coupling parameters in the generalized Wigner sur- 
mise and in the large mixed matrices are related via the 
local eigenvalue density of the latter. This relation is well 
approximated by Eq. ( |3.7[ ). 

We expect that the results derived in this paper will 
be useful in numerical and/or experimental studies of 
systems with mixed symmetries. 
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Appendix A: Analysis of the Gibbs Phenomenon 

The spacing distribution of ensembles interpolating be- 
tween Poisson and RMT reveal a Gibbs-like phenomenon 
close to the Poisson limit, i.e., for small A: P{s;X) does 
not converge uniformly to the Poisson curve e~* at s = 0. 
Rather, there is an overshoot whose amount does not van- 
ish in the A — >■ limit and whose position s approaches 
in this limit. In this appendix we work out the value 
and position of this maximum. We start with a brief re- 
view of the Gibbs phenomenon in the Fourier transform, 
as known from textbooks such as [26] (which, however, 
mostly discuss the Gibbs phenomenon only in the Fourier 
series) . 



1. Gibbs Phenomenon in the Fourier transform 

The Gibbs phenomenon is related to the convergence of 
the inverse Fourier transform with a cut-off in the integral 
(or to the convergence of the Fourier series with a cut- 
off in the sum) towards the original function /. Let us 
denote its Fourier transform by F, 



of the original function with the Dirichlet kernel 



F{oj) 



1 

2n 



(Al) 



and the result of the cut-off inverse transform by / with 
two arguments 

/(s;A) = / duje"^'F{uj). (A2) 
J-i/x 

These formulas can be combined into a convolution 

/>oo 

/(s;A)=/ ds'f{s')Sx{s-s') (A3) 



(5a (s - s') = 



2tt 



l/A 



d^^Ms-s') ^ sin[(s - s')/X 



-i/x 



tt{s - 



(A4) 



The question is how in the A — )• limii[^ /(s; A) is re- 
lated to the original function f{s). If /(s) is smooth and 
absolutely integrable, /(s; A) approaches it everywhere. 
Accordingly, the Dirichlet kernel approaches the delta 
distribution in the sense of acting on smooth test func- 
tions. 

At discontinuities of the original function /(s), how- 
ever, /(s; A) approaches the average of the left and right 
limit of f{s). Intuitively, this comes from the nonzero 
width of the Dirichlet kernel]^ which in the convolution 
(A3 ) probes both sides of the discontinuity. Furthermore, 



the functions /(s; A) for fixed A possess maxima and min- 
ima whose positions move, in the limit A — >■ 0, towards 
the discontinuity and whose values over- and undershoot 
the function. This is the Gibbs phenomenon. 

For dcfinitcncss let us consider a set of exponentially 
decaying functions with a jump discontinuity of unit sizaj 
at s = 0, 



/(«) 



for s < , 
e-'*" for s > , 



(A5) 



that include the Poisson curve (d = 1) and the Heaviside 
function (d = 0) . The Fourier transforms are 



Fiuj) 



1 



1 



2tt d + iuj 

and the cut-off inverse transforms read 



(A6) 



fis;X) 



2n 



Ei (ds - ij^ - Ei (ds + ij^ 



(A7) 

For the Poisson curve these functions are plotted for three 
small values of A in Fig. [12] (top), where several maxima 
above and several minima below are clearly visible. 

To analyze the limit A — > we can zoom into the region 
of small s, of size proportional to A. This amounts to 
considering functions of a rescaled argument 



s 



(A8) 



in a constant s-range. We define 
/(S;A) = /(.?A;A) 



^ The A — > limit of f{s; A) is also denoted as the principal value. 
* This nonzero width is relevant in many areas of physics such as 

band-limited signals, ringing, and diffraction of waves at slits. 
® As all formulas are linear in f(s), the case of arbitrary jumps is 

completely analogous. 
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FIG. 12. The Gibbs phenomenon in the Fourier transform of 
the Poisson curve, with A — 0.01,0.005,0.0025, respectively. 
Top: the functions /(s; A) approaching the original function 
e~" (dashed), first maximum moving to the left as A decreases. 
Bottom: the rescaled functions f{s; A) approaching the limit- 
ing function g{s) (dashed) with decreasing A, see text. 



-dsX 



27r 



[Ei(dSA 



Ei{dSX + iS)]. (A9) 



Keeping s fixed, these functions have a well-defined limit 
A ^ 0, 



g{S) = lini /(S; A) 



[Ei(-zS) - Ei(zS)] = I 



Si(S) 



(AlO) 

with the sine integral Si(s) = l^dx smxjx. As Fig. 12 
(bottom) shows, this limiting function captures infinitely 
many maxima at s = tt, Stt, . . . and infinitely many min- 
ima at s = 27r, 47r, . . . The overshoot at the first maximum 
is the well-known number 



Si(7r) 



1 = 0.0894899. 



(All) 



Concerning the convergence of the Fourier transform, we 
conclude that in the limit A — >■ the functions /(s; A) 
have a maximum at s = ttA, with an overshoot approach- 
ing 8.9%. 

Note that the limiting function g is the same for all 
these functions independently of the decay constants d, 
i.e., it is solely determined by the discontinuity. In other 
words, the smooth part of the function /(s) drops out 
when going from /(s; A) to g(s) in the A — > limit, see 
Eq. ([^ vs ( |ATo| . 

This can be shown to be universal. Rescaling the inte- 



gration variable in Eq. (A3 1 and using \6\(\x) — 8\(x) 
one has 

f{s; A) = ds"f{Xs")6iq - s") , (A12) 

/(S;A)= / ds"/(As")'5i(S-s"), (A13) 
Jo 

where we still assume f{s < 0) = for simplicity. The 
limiting function is 



ff(S) = /(0+)/ diJi(t) = /(0+) 



Si(5) 



(A14) 



which agrees with (AlO) for all functions with /(O''") = 1. 



The last equation in particular relates the Dirichlet 
kernel Sx and the limiting function g. Therefore, /(s; A) 
can also be reconstructed by a convolution with (the 
derivative of) g, 



/(s;A) 



ds'fis') 



(A15) 



where g'{x) = dg/dx. 



2. Gibbs Phenomenon in the Poisson to RMT 
transitions 

For the Gibbs-like phenomenon in the mixed spacing 
distributions we start with the Poisson to GSE case. 
For this transition we found the spacing distribution 
Eq. (2.26). In analogy to the previous subsection we 



rescale the argument and define 



Po^4(S; A) = Po^4(SA; A) = CX^S^e 

2\ {xD\s)^+2X^xDs 



X 



dx{l-x^)e^ 



(A16) 
erfc(a::_DAs + A) . 



-1 



In the limit A — > we make use of the behavior of C(A) 
and ^(A), 



D{X)^^ and C(A)~^, (A17) 



to arrive at the limiting function 
go^4{s) = hm Po^iis; A) 

A — ^0 
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This concludes our empirical results on the Gibbs-like 
phenomenon. 

Concerning the analogies at a more fundamental level, 
the spacing distribution Po^i{s; A) is related to the inte- 
gral (^231) 
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FIG. 13. The rescaled spacing distributions Po_s-4(s; A) 
in the Poisson to GSE transition, Eq. (A16I, for A = 



0.05, 0.025, 0.01 (maxima increasing) approaching the limit- 
ing function go-^4(s), Eq. (A18I (dashed). 



jI{S/X)=J^ dpe-PSx{S,p) (A21) 

of the unperturbed Poisson distribution with the kernel 

SxiS,p) = J da...dc3Pa{a)... PcM) (A22) 

x6(^S- \^(a-b-p/Xy+4{c^c^)^ . 

The nonzero width of this kernel causes the Gibbs phe- 
nomenon in the spacing distribution near the disconti- 
nuity of the Poisson distribution at p = 0. Note 
that in the limit A — ?> the second line of Eq. (A22) 
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s 



approaches S{S — p), thus decoupling from the integrals 
over a, . . . , C3. The latter are normalized by construction 
so that the kernel 6\{S,p) approaches 6{S — p). 

There are (at least) two features that are different from 
the Fourier case. First, the kernel is not a function of S — 



p, and thus Eq. ( A21 ) is not a convolution, in contrast to 



r, , ~2\ /— -sl nf~ir,\ r,~] f \ T Q\ thc Fourlcr casc, Eq. ( Al 5 1 . Second, at the discontinuity 
Z + s ),/T:e evti[s/Z) - Is . (AiS) 



In Fig. 13 we plot this function, together with Po--i-4(s; A) 



approaching it for small A. 

Again, this limiting function captures a maximum, 
which can numerically be determined to be at s — 
3.76023 with a value of 1.43453. As before, the phe- 
nomenon is solely determined by the discontinuity of the 
Poisson curve at s = as the original Poisson distribu- 



^b-i-4(0; A) = is not the average (equal to 1/2) of the 
left and right limit of the original Poisson curve e~'^ (put 
to zero for negative p). 



Appendix B: Explicit calculation of spacing 
distributions 



Poisson to GSE 



tion e ^ can be shown to drop out from Eqs. ( A16| ) to 



(A18) 



So in the transition from Poisson to GSE, Po->4(s;A) 
in the limit of small A has a maximum at s = 3.76A 
overshooting the Poissonian e^" by 43.5%. Likewise, in 
the other transitions Pq^i{s\X) and Po-!-2(s; A) we have 
equivalently defined limiting functions 



5o^i(-?) = ^.5e-^"/8/o(SV8) , (A19) 



We start from Eq. ( 2.25 1, transform co, . . . , C3 to polar 



50^2(5) 



erfi(5/2) 



(A20) 



These have maxima at s = 2.51 and s = 3.00, overshoot- 
ing the Poisson curve by 17.5% and 28.5%, respectively, 
as quoted in the body of the paper. We observe that these 
numbers grow with the Dyson index /?' of the perturbing 
ensemble. 



coordinates with — c^c^, and introduce u = a + b and 
t = a — b. This yields 

I{S)oc dpdcc' dwdte-P-^("'+*')+^-^-^' 

JQ J -00 

xS(s- ^/t^ + 4c2) 

"■^ ' die ^ 2A I 

i 

s(S-Vt^+AcA, (Bl) 



oc / dpdc c 

'0 J -00 



where in the last step we have integrated over u. We now 
use the (5-function to integrate over c, resulting in 



From the small-s behavior go^4{s) = s^/12- 



-0(s6^ 



we 



conclude that Po^4(s; A) = s*/{l2\'^) + 0{s^) for small 



I{S) oc / dt (S^ - St^) e 



dp e 



p-'+2Xp(2X + t) 



A, which reproduces our observation in Eqs. (2.301 and 



oc S*^ e 



(2.31). Analogous agreement is obtained with Eqs. (2.91, 



(2.10) and Eqs. (2.19), (2.201 for the other two cases 





2,2 



dx {I — X )e 1- 



AS): 



(B2) 
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where we have integrated over p and substituted t = xS. 

T he p roper level spacing distribution Po-j-4(5), 
Eq. (2.26), is obtained from J{S) by rescaling and nor- 
malization, i.e., Po^4(s) = C ■ J{2Ds)/{2D)'^. Defining 
the moments of the distribution, 



T„= / dS'5"J(S'), 
/o 



we obtain from Eq. (1.2) 



D^^ and C=(^ 
Explicit evaluation of Iq and Ii gives 

16 



T 16 ^ 



Ii = 







(4a;3 + 2.T)e-^ + ^(4a;4 +4a;2 - 1) erf (x) 



(B3) 



(B4) 



(B5) 



(B6) 



from which we obtain Eqs. (2.27) and (2.28). 



GOE to GSE 



We consider the matrix H in Eq. (2.40). With a small 
change in notation for Hi, we have 



H ^ Hi (E)l2 + XH4 



/A C 0' 
A C 
C B 

,0 C By 



+ A 



a Co + ci + ic2^ 

a — ci + ic2 cq — ic^ 

Co — ic3 —ci ~ ic2 b 

\ci — ic2 Co + ic3 b 



(B7) 



where the variances of the random variables are given 
by Eq. (2.2). If two variables are Gaussian distributed 
with variances and cr^, their sum is again Gaussian 
distributed with variance + cr|. Since H depends on 
A,B,C and a, 6, co only through the combinations A + 
Xa, B + Xb, C + Acq, we can immediately integrate out 
A, B, C, with the corresponding change in the variances 
of a, 6, Cq. To simplify the notation, we absorb A in H4 
and divide all matrix elements of H b y \/l + A^. This 
yields a problem equivalent to Eq. (B7), 



H 



with 



/a Co + ic3 ci + JC2> 

a —ci + ic2 Co — ic3 

Co ~ ic3 — Ci — ic2 b 

ci - ic2 Co + ic3 b 



(B8) 



2a.^. - 



A2 



1 + A2 



(B9) 



where i = 1,2,3. The matrix in Eq. (B8) has two non- 



degenerate eigenvalues whose spacing is given by 



S = 



1/2 



(BIO) 



where we have again written S instead of s since we still 
need to enforce the normalizations (1.2). The spacing 



distribution is proportional to the integral 



liS) 



da db dcf)dcidc2dcz e 



(Bll) 



where repeated indices indicate a sum over i and j from 1 
to 3. We now transform ci, C2, C3 to spherical coordinates 
with = CiCi and introduce u = a + b and t = a — b. 
This yields 



7(5) oc / dT/didcodcc^e- J +4^0)-^ 



X 5 ( S- ^|t^^+Acl + Ac^^ ) . 



(B12) 



The integral over u can be performed trivially and only 
results in a prefactor. Using the (5-function to integrate 
over c, we obtain 

/>oo 

7(5) oc / dtdcoe-j(*'+4^«)-l^(^'-*'-4^o) 



X - t2 - 4c2 e {S^ -t^- 4cl) , (B13) 

where we have used the symmetries of the integrand to 
raise the lower limit of the integrations to zero. We now 
perform the transformation 



t = Sx and co = -Syy/Y 



(B14) 



with Jacobian — x'^ . Since t and co are non- 



negative, so are x and y. The 0-function in Eq. (B13) 



then implies < a;,j/ < 1. Reinserting the definition of 
from Eq. (B9), we obtain 



I{S)(xS^ I dxdy{\-x^)y^l-y'^ 

,[A^ + (l-:r^)(l-^S,2)] 



X e -iA-! 



(B15) 



We now substitute y = cos note that cos^ (j) = 
|(1 — cos 20), and use the integral representation [18l 
Eq. (9.6.19)] of the modified Bessel functions Iq and Ii 
to obtain after some algebra 



I{S) oc S'^e-^^ 



dx{l~x^)ei^ [Io{z) - h{z)] 

(B16) 
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with z = (1 — x'^)S'^ /{SX"^). This corresponds to 
Eq. ( |2.41| with S = v^ADs. The properly normal- 
ized spacing distribution is therefore given by Pi^i{s) = 
CJ{V8XDs)/{V8XD)^. Defining 



T = 



dSS^'JiS) 



we obtain from Eq. (1.2) 



D 



and C 



iVsxDY^ 



Explicit evaluation of Xq and Ii gives 
A2 



Ii = 16A^ 



1 + A^ 
A(l-A^) 

(1 + A2)2 



arccot A 



(B17) 

(B18) 

(B19) 
(B20) 



from which we obtain Eqs. (2.42) and (2.43) 



Appendix C: Perturbation of a large GSE matrix by 
a non-self-dual matrix 



We consider a mixed 2N x 2A'' matrix that interpo- 
lates between the GSE and one of the other Gaussian 
ensembles, 



H - Hi^ J-— Hp, 



(CI) 



where i?4 is taken from the GSE and Hpi from the GOE 
or GUE. We study this matrix for large N in first-order 
degenerate perturbation theory to show similarities be- 
tween the two different perturbations and to make a con- 
nection to the case of GSE to non-self-dual GUE for 
N = 2, which was treated in Sec. |II H[ 

Degenerate perturbation theory predicts that each of 
the N previously degenerate eigenvalue pairs splits up 
and that the shifts of the two members of the pair are 
the eigenvalues of the matrix 



A 

P4(0)S4 



My with My = (V^, \Hp, I V,> ; i, J = 1, 2 . (C2) 



The |'0i,2) are the orthonormal eigenvectors of the un- 
perturbed matrix that span the degenerate subspace 
of the eigenvalue pair under consideration. 

We show in the following that M is a 2 x 2 GUE matrix 
for /?' = 2 as well as for = 1, in the latter case with a 



normalization different from Eq. (2.2). 



1. GUE 

This case is very simple, because the GUE is invariant 
under unitary transformations, which contain the sym- 
plectic transformations. This means that the transforma- 
tion diagonalizing the GSE matrix can be absorbed 



in H2 without loss of generality, and therefore one can 
choose \ipi)k = 5ik with i = 1, 2 and fc = 1, . . . , 27V. Thus, 
we obtain 



2N 



M,, = {H2)ki hj = {H2) 



(C3) 



which is obviously a 2 x 2 matrix from the GUE with 
the usual normalization, Eq. (2.2). As this holds also 
for iV = 2, it is a perturbative explanation for the fact 
that in the limit A — > the spacings between previously 
degenerate eigenvalues are distributed exactly like the 
ones of 2 x 2 GUE matrices. 



2. GOE 

We will show that in this case M is again a matrix from 
the GUE with the only difference that the variances of 
its elements are only half as large as in the previous sub- 
section. This case is a bit more involved because one 
cannot generally diagonalize a self-dual matrix by an or- 
thogonal transformation (which would preserve the prob- 
ability distribution of Hi), and thus it is impossible to 
choose the eigenvectors of H4 as in the previous subsec- 
tion. Explicitly, the matrix elements read 



My- = {^p^\Hl\lPJ 



+z((^rii?ii^f)-(V'ri^iiV'r)), 

where we split the eigenvectors l^pi) in real and imaginary 
parts: = \ipl^) + i\ipi^), and Hi is real. 

We will now show that the four vectors IV'™)) IV'i™); 
and IV'™) are orthogonal in the limit of infinite 
matrix size. For some combinations of them one can show 
this also for finite N using the quaternionic structure of 
the eigenvectors. 



(V'lh _ 
(^2 



9l 92 



In) , 



with quaternions in matrix representation 

/ (0) , ■ (3) (1) , . (2)^ 
Ik - \ (1) , . (2) (0) . (3) 



(C5) 



(C6) 



One can read off immediately that 

(V'riV'2°) = (V'ri^n = o, (c?) 

= (C8) 
(V'riV'r) = (^riV'2^) , (C9) 

i.e., there are only two independent scalar products. 

Let us assume that for large TV the q^^^ can be treated 
as independent random variables with mean value zero. 
Then the mean values of those scalar products are zero 
as well, e.g., 

N 

((V'ri^r)) = E('^i^r^+«i^i'^) = 0' (CIO) 
fc=i 
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where the outer angular brackets indicate an average over 
the random matrix ensemble. From the normalization of 
the eigenvectors \'ipi) the variances of the g^''^ are propor- 
tional to This yields for the variances of the scalar 
products 



sembles H, 



N 

E 

k=l 

N 



)) 



E- 

^ N 



k=l 



1 

N 



(Cll) 



and likewise for (tA^^"!?/;™). Since in the N oo limit 
both the mean values and the variances of the scalar 
products vanish, the four vectors become orthogonal in 
this limit for every single realization of the random ma- 
trix. We have checked this numerically, which implies 
that the assumption of the independence of the g^''^ was 
valid. 

As for the normalization of the four vectors, the 
squared norms of the real and imaginary parts agree on 
average and sum up to 1 due to the normalization of the 
eigenvectors Invoking the central limit theorem, we 
observe that in the limit N ^ oo the norms of the real 
and imaginary parts equal 1 / V2 even for a single realiza- 
tion of the random matrix. Hence, multiplying the four 
real vectors Itpl"), li/'l™), |?A2°)> and \tp'^) by ^/2, one ob- 
tains, in the limit N oo, an orthonormal real basis in 
the subspace under consideration. 

Finally, we use the fact that the matrix elements of 
a GOE matrix Hi are independent random numbers in 
every orthonormal (real) basis, with variances 1 and 1/2 
on and off the diagonal, respectively. Thus we conclude 
that the M^- are also independent random numbers with 
variances 

([Af-f ) = ([M-]^) = 2, (C12) 
([M-f ) = ([Afl-]') = i. (C13) 

These are half the variances of a GUE matrix, which 
is equivalent to a multiplication of each element of M 
by l/\/2- This explains the rcscaling of the coupling pa- 
rameter in the definitions of Pl^i{si; A) and Pi^i{s2', A), 



Eqs. ( [3l8l ) and ( |3l9 |. 

For small TV the argumentation in this section does 
not work. Presumably, this is the reason why the spac- 
ing distributions for the transition from GSE to GOE 
differ from those for the transition from GSE to GUE in 
the case of 4 x 4 matrices (not shown in this paper, but 
checked numerically), whereas they match very well for 
large matrices. 



Appendix D: Perturbative calculation of the relation 
between eigenvalue density and coupling parameter 

We consider a diagonal Poissonian matrix Hq per- 
turbed by a matrix taken from one of the Gaussian en- 



H = Ha + aH, 



(Dl) 



where Hpi is chosen in the usual normalization, see 
Eq. (2.2). The calculations are done for arbitrary matrix 
dimension, which will be sent to infinity at the end. We 
denote the number of generically non-degenerate eigen- 
values by N , i.e., we consider N x N matrices. If Hpi 
is taken from the GSE, these are quaternion valued and 
correspond to complex 2A^ x 27V matrices. 

To obtain an iV-independent eigenvalue density of the 
Poissonian ensemble, we define the probability distribu- 
tion of the individual eigenvalues 9i of Hq by 



(D2) 



where Vq is some iV-independent probability distribu- 
tion. Both Vq and Vq are normalized to one. The eigen- 
value density of the Poissonian ensemble is thus 

Po{9) = NVoio) = Me/N) = Me) , (d3) 

where we have defined 9 ~ 9/N. Generically we have 
9, = 0{N) and 9, = 0(1). 

We now consider a fixed spacing S between two adja- 
cent eigenvalues of iJo, and ^2=^*1+ S. The remain- 
ing eigenvalues have to reside outside the interval {9i, 92). 
This results in the conditional probability distribution 



V°Q 



1 



)^Ti'Pr{o,iN) 



N 







•Po(9) 



for 9, e (01, 
otherwise . 



(D4) 



The eigenvalue density is assumed to be almost unaf- 
fected by the perturbation, which is confirmed in Fig. [6] 
(top). Of course, this assumption is expected to hold 
only for small values of the coupling parameter. 

We want to calculate the effect of the perturbation 
on the spacing S. If the remaining eigenvalues of Hq 
are close to 9i or 02 we have to apply almost-degenerate 
perturbation theory. Up to second order in a we obtain 
for the perturbation of the spacing 



AS* = (^EVD 



(Ho + aHp,) 



first-order almost-dcgcncratc pert, theory 
, 2 



N 



a 



i=3\ei<^w 



(D5) 



second-order perturbation theory 



where the absolute values are taken with respect to 
the real/complex/quaternionic standard norm, EVD de- 
notes the difference of the two eigenvalues of the matrix 
{Hq -j- aHpi)i^i^ that correspond to the unperturbed eigen- 
values, and W is the interval in which eigenvalues have 
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to be considered almost degenerate with 9i or 02- This 
is defined by the eigenvalue range {9i — Cw,&2 + Cw), 
where we choose Cw = cj^^iV^a with < e < 1 and 

> I. This choice ensures that the closest possible 
eigenvalue outside W cannot give a second-order contri- 
bution of lower order in a than the almost-degenerate 
partj^ Note that the "degenerate window" W grows 
with N. Therefore arbitrarily distant eigenvalues are con- 
sidered almost degenerate in the limit N ^ oo, which is 
justified because almost-degenerate perturbation theory 
is valid for any difference of eigenvalues. 

Considering the first-order contribution, we have to 
deal with the matrix 

Mki = (i/o)« + a{Hp')ki = OkSki + a{Hp.)ki , (D6) 

where the indices k and I run over all values for which 
the eigenvalues Ok and 9i are localized in W, which in- 
cludes at least 0i and 62 ■ This is a matrix taken from the 
Poissonian ensemble perturbed by a matrix taken from 
one of the Gaussian ensembles, but unlike H defined in 



Eq. (Dl I it has a constant eigenvalue density in the limit 



N — T' 00. To show this, we first consider the density at 
the lower end of the interval W, 



lim NVoiOi-Cw) 

N-~>-oo 



hm A i ^ 



lim Vo 6*1 



N 



Cw 



(D7) 



Pol 



This is the same as the eigenvalue density at the other 
end of W, 



lim NVo{02 + Cw) 



lim Vq 6*1 



s_ 

N 



lim Vo 



N 



S + C 



w 



N 

VoiOi) 



(D8) 



PoiOi). 



Thus the spectrum of M can be unfolded by multiplying 
with the local eigenvalue density. 



Po{Oi) Mki ^ PoiOiWkSm 

unfolded 



+ apo{ei)iHp,)ki. (D9) 

effective coupling 



Therefore we can define a new effective coupling param- 
eter that solely determines the magnitude of the pertur- 
bation as in Sec. IIII Al 



The second-order contribution to AS* in Eq. (D5) is a 



sum of at most TV — 2 independent random numbers. As 
all of these random numbers have the same distribution 
we pick out one of them, 



a 



with e^^W, (DIO) 



where we defined a = |(_ff^/)ii| and b 
probability distribution is given by 



Its 



61—Civ 



:+Cw 



1 



d9j^Vrie/N) (Dll) 



X / dadhVi3'{a)Vfi'{b)5 



X — a 



h-e Oi- 



where we renamed Oi — 9 for convenience. The distribu- 
tion Vi3' depends on the symmetry class of the perturbing 
ensemble (a and b are squared sums of f3' Gaussian ran- 
dom variables) . The moments of this distribution are 



After a short calculation, we obtain 



Pm = 



(D12) 




(D13) 



In the limit iV — > 00, all terms that are divided by N 
in the arguments of 'P™* can be neglected. This can be 
done in spite of being integrated to 00, because the 
last part of the integrand (in square brackets) suppresses 
the large-6' region and because Pq"* is a probability den- 
sity that has to converge to for large argument. Also, 
limAr^oo^r*(^i) = ^o(^i) = Po{di)- We thus obtain 



2po{9i) 
N 



d9 / dadbVf}'{a)Vf}>{b) (D14) 
-00 Jo 



a 



Let us denote the second line of Eq. (D5 1 by AS**^^-*. It is 



0{Npi), and therefore its mean value becomes zero for 
iV — 00, as it is suppressed by . The same holds 
for the second moment of A5*^^-', which goes like N^^' . 
Thus the distribution of AS'^^^ is a delta function at zero, 
and we can neglect its contribution to the perturbation 
of the spacing. The linear relation between eigenvalue 
density and coupling parameter could hence be shown 
up to second-order perturbation theory. 



Appendix E: Method for fits to the surmises 

An eigenvalue 9i = 62 + Cw at the border of W yields 

a second-order shift of A02 = -a^\{H^,),2\V{C^^^ N^a) = gj^^^g ^^^^ analytical formulas for the small ma- 

-a\{Hisi)i2\ /(C^yN'^). trices contain integrals, it takes some time to compute 
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them numerically. In order to get good fits to data in 
a reasonable time, a list of 1000 A-values in the interval 
(0.01, 10) was created, with 



A, = 0.01 • 1000959 ; i = l, 



1000 . 



(El) 



For each and each surmise, the corresponding spacing 
distribution was stored. The pure cases A = and A = oo 
were included as well. 

As a measure of the fit quality, we use the L2-distance 



dx [f{x) - g(x)]^ 



1/2 



(E2) 



between the fit and the numerical data. The fitting was 
done by calculating the A2 value of each spacing distri- 
bution in the list. From the one resulting in the smallest 
A2 we read off the coupling A. Note that the largest A2 
we encounter in all the fits is 0.019. For comparison, the 
L2-norms of the pure Wigner surmises Pp{s) range from 
0.71 for /3 = to 0.94 for /3 = 4. We give no error bars, 
because the statistical errors of A obtained by methods 
such as Jackknife were negligibly small. This is also the 
reason why we use A2 instead of a statistical quantity 
like chi-squared as a measure of the fit quality. 



Appendix F: Construction of a self-dual GUE 

In the following we construct a Hermitian, self-dual 
2N X 27V matrix whose eigenvalues are twofold degener- 
ate and whose non-degenerate eigenvalues correspond to 
those of a matrix from the GUE. We start with a ma- 
trix M that contains an x iV GUE matrix H and its 
complex conjugate (equal to the transpose). 



M = 



H Oat 



(Fl) 



The eigenvalues of M are obviously those of H, but now 
twofold degenerate as desired. However, M is not self- 
dual. To transform M into a self-dual matrix without 
changing its eigenvalues, we apply an orthogonal trans- 
formation 



O = 



/ 1 
/ 

0010 
0000 



0000 
0000 

0100 
0000 
0001 
0000 



\ 
\o 










1 









1 




0000 
1000 
0000 
0010 



0000 
0000 

0000 
0100 
0000 
0001 



0000 
0000 



0\ 

\ 









1 








0/ 

01/ 



= 0' =0- 



(F2) 

which transforms a matrix by exchanging every 27i-th 
row and column with the (A'' + 2n — l)-th one. Each of 
the four blocks is a square matrix of dimension A^. This 



is in complete analogy to the construction of a self-dual 
4x4 GUE matrix in Sec. HTGI 

We now show that the transformed matrix MO is 
self-dual, the condition for which is 

M = OJOM^OJ^O (F3) 



with 



J = l 



N ' 



-1 

1 



(F4) 



Multiplying J by O from the left and the right in- 
terchanges the second, forth, . . . with the (A^ -I- l)-th, 
(A^ -I- 3)-th, . . . column and row. We thus obtain 



f&N -In 



OJO = 

= -OJO = 



Oat Ijv 



(F5) 



and thus 

OJOM^OJ'^0 = 









Oat 


H 









= M, 



H 



(F6) 



which proves Eq. p^. O'^MO can therefore be writ- 
ten as a quaternion matrix with real quaternions and 
their conjugates at the transposed position. Each of these 
quaternions stands for a matrix of the form 



Co + ic3 ci + ic2 
-Ci+ ic2 Co — ic3 



q p 

-p* q* 



(F7) 



with complex numbers p and q. Evidently, p has to be 
zero for each quaternion in O^ MO, because our original 
M generically contains no element which is the negative 
complex conjugate of any other, and we only exchanged 
elements by applying O. This means that at least half 
of the matrix elements are zero. In the original M, ex- 
actly half of the matrix elements were zero, while the 
other half were random variables which depended on a 
total of A^^ real parameters, so the same has to hold for 
O^ MO. From this and Hermiticity if follows that every 
off-diagonal q has to be an independent complex random 
number, while the q on the diagonal are real, so that 
there are again A^^ real degrees of freedom. 

With this equivalence proven, one can construct a self- 
dual GUE matrix by taking a matrix from the GSE and 
set its ci and C2 components to zero. This matrix has 
the same joint probability density of the eigenvalues as 
a.n NxN matrix taken from the GUE, as it is related to a 
matrix of the form of M by a fixed basis transformation. 
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